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Abstract—The evolution of space missions and related systems has been 
promoting innovation and creation of ideas for new technologies for 
decades. The search for innovative solutions that combine the optimization of 
resources and materials guide the recent research in the space area. The 
objective of this study is to analyze the behavior of two bodies connected by 
tethers and can be a solution for reducing costs in space missions, one of the 
concepts that have the potential to fulfill the objective of efficient transport 
space. In this paper it is discussed the motion of two massive bodies 
connected by tethers in keplerian motion in a central force field, their 
viability, the rotational dynamics and the system behavior in the space 
environment. Models will be created that simulate and explain the dynamics 
of the object and that analyzes the main parameters for determining the 


mathematical modeling. 





stability and the uniform rotations conditions. 





I. INTRODUCTION 


Tethered systems have many areas of application and has 
been studied in several published articles ({1]-[9]). 
Beletsky and Levin[1], begins by setting the scene for 
tethers in space summarizing possible applications and 
also discussing fact and fiction, analyzing clearly the main 
parameters and applications for Tethers Systems, as the 
density of the material the effective forces, orbital 
dynamics, mechanics models, attitude and possible 
disturbances for a flexible tethers with end masses, 


massless and massive variations. 


The motion of tethers considering dumbbell 
oscillations, bodies in the central field of vibrating forces 
(Burov et al. ([10]-[14]), showing stability solutions for 
angles on the chaotic dynamics in elliptical orbit, 
analyzing the problem in another aspect of Moon-tethered 
pendulum, to considering the uniform rotations of a two- 
body tethered system in planar motion and the control the 
length of the tether [15]. In studies ({15]-[21]) were 
suggested methods of controlling the geometric 
configurationand Moon-tethered system with variable 


tether length in restricted three-body problem ([16] and 
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[22]) were studied and demonstrated important results for 
tethers applications in systems. 


II. DYNAMICS OF THE PROBLEM 


Consider planar motion of a dumbbell body in the central 
Newtonian gravitational field (Fig.1). The tether length (/) 
is small compared to the orbit the system center of mass 
(cmg), moves along an elliptic keplerian orbit. 


y 4 











Fig.l: Geometry of the tether system in reference frame. 
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The body consists of two-point masses (m, and m3) 
connected by a light tether 


p 


~ T+e cos(v) (1) 


p 


where p is the focal parameter, e the eccentricity and v the 
true anomaly of the orbit. The system coordinates are 
Xo = pcos(v) 
Yo = psin(v) 

Xı = Xo + L cos(v + o) 

Yı = Yo +1, sin + o) 

X> = Xo — l, cos(v + ø) 

Y2 = Yo — hz sin + @) 


(2) 


1.1. Potential and Kinetic Energy 
The potential energy of the system can be obtained by the 


following expression, 
Hoa _ Honn 
ml l 








V=- (3) 


where m; represents the mass of the point; 7;the position of 
the point mass with respect to center of the Earth; Uy = 
GM; G is the universal gravitational constant and M the 
mass of the Earth. Introduced the parameters, u and m, 


given by 
p= (4) 
m=m +m, (5) 
Substitute on Eq. (3)leads to 
V = -uoy (—" + =) 6) 
vxi ty x2 +y 


The Eq. (7)is cumbersome and can be simplified, 
introducing a new parameter A = = and assuming A « 1 
the tether length is much smaller than the focal parameter 


p. The Taylor Series expansion up to the 2” order of A for 
the potential energy is 


V 
_ MUp(1 + e cos(v)) 
g p (7) 
F m(u — 1)uuo(1 + e cos(v))? (1 + 3 cos(2@))A? 
4p 


The kinetic energy of the system can be written as 
_ mvi? ~ m,|v2\° 


ee 8 
T . 5 (8) 
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"N 1 [2p?(1 + e? + 2e cos(v))v? 
~ 4 (1 + ecos(v))* 


—2u-De(?PW+¢)? 9) 


+ 12) 


1.2. Lagrange Equations of Motion 
For the subsequent analysis, the generalized coordinates @ 


and l are used and the system is assumed to be subject to 
the gravity-gradient forces. 


d (ðL) dL 
Tg) dp (10) 
d ðL dL 
aa 
3 oes at 
l (Suosin(29)(1 + ecos(v)) + 2p? + @)) (11) 


+ 4p3l(v + @) = 0 
The equation is rewritten as a function of the true 
anomaly v (Eq. (11)- (12)) 


oo 12 
( ) T; (12) 
d d d 
sees eee 2 — 13 
T Wo (1 + e cos(v)) p (13) 
with 
mo (14) 
0 p3 


The equations of motion of the spacecraft can be 
written as 


(1 + ecos(v))o” 


+2 G (1 + ecos(v)) 


(15) 
—e sin) (o' +1) 
+ 3cos(@) sin(g) = 0 
considering a movement with uniform rotations, 
p = wV + Po (16) 
l(v) = nv) —— 17 
ay) 1 + ecos(v) oe 


Suppose the following relation for the tether 
performance (Eq. \ref{equal6}), the Eq. \ref{equal4} with 
respect to the true anomaly v takes the form: 


n' (v) B 3 sin(2(wv + Po)) 
nv) 4(w +1)(1 + e cos(v)) 


The analytical integration function is difficult but 


(18) 


substituting values for the variable w and solving the 
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equation with respect to the variable 7, it is possible to 
obtain the following solutions (Table 1 and Figures 
\ref{fig2} - \ref{fig5}). Analytical solutions were found 
for this w values previously chosen. The values for w = 


(+2, E, $a: pL pL pa +=, +>) closed-form no 


solution. For w = 0 was examined in previous studies [12] 
where the relative equilibrium was studied. 


@w=1 


= 0.04 
= 0.2 
= 0.3 
=0.5 
= 0.8 


o oo 8 ® 








Fig.2: Tether length control for w = 1. 
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Fig.3: Tether length control for w = 2. 


a= -4 
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Figure 4 -Tether length control for w = —4 


The control laws are periodic in true anomaly (v). In 
Fig.2to Figure 4describe the change in tether length 
depending for different eccentricities,e = 
(0.04, 0.2, 0.3, 0.5, 0.8 ). The formulation also guarantees 
uniform rotation for fractional values of w, as seen in 
Figures6—9. 
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Fig.5: Tether length control for w = -+ 
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Fig.6: Tether length control for w = > 


1.3. Uniform Rotations: @ = wv + Po 

Substitute (p = Wv + Qo) to Eq. (14)) and considering 
Eq. (17), it is possible to obtain the nonlinear differential 
equation that describes the motion of the tether. 


A general analytical integration for all values of w 
cannot be found but for some values w the closed form 
solutions obtain. A numerical integration has been 
performed, where w and @p, were substituted previously, 


with the objective of analyzing the behavior of the cable 
system. The logarithmic plots of In (=) for a number of 


0 
fractional angular velocities (Figures 6 - 9) was made. For 


the particular case when the value of w is zero (w = 0) 


and pọ is constant, there are a uniform for o = — k= 
{0,1,2,3,4,...}. 
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Fig.7: The control law for the tether length for w = 0 and 
Po = 0. 


An approximate solution for small eccentricities can be 
obtained using Taylor series in Eq. (17)and obtain 
applying the series of order 3, in the variable e obtain: 

n'(v) 3e? cos?(v) sin(2vw) 

n) 4@ +) 

3e? cos? (v) sin(2vw) 
g 4(w + 1) 
3e cos(v) sin(2vw) 
4(w + 1) 
3 sin(2vw) 
= A(o + 1) 





(19) 


One obtains the expression that has analytical 


integration, however’ this procedure introduces 


impossibilities in the system and restricts the solutions to 
the variable $\omega$, which after integration generates 
singularities for values: w= (- =, —1,- =, 0,-, 1, =), 
Increasing the expansion terms for order 5, the 
also increase, w = 


5 3 1 1 3 5 
(-=, =A; ag = 1, g) 0,-, 1,>,2,2). 
2 2 2 2 2 2 


singularities 


1.4. Stability Conditions 
We study stability using to the Floquet Theory, let’s 


linearize the equation of motion in the vicinity of solution 
(Eq. (16) ), it is possible to check the stability analyzing 
system behavior around small oscillations. Considering 
Eq. (19)andEq. (20)analyzing the neighborhood of the 
point Qo = 0, it follows: 


p = wv + do (20) 
Substituting Eq. (17)in Eq.(15)one obtains: 
; n'v) 
4(w + 1+8) + e cos(v)) —— 
nv) 


+3 sin(2(wv + 5@)) a 


+ 2(1 + ecos(v))6@”" = 0 
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Replacing Eq. (18), concerning the n(v) variable, 
d''(1 + ecos(v)) 

+ 3cos(2wv + dq) sin(dq@) 

3sin(6~") _ 

2(w+1) 


(22) 
0 


This is the nonlinear equation of the perturbed motion 
\cite{ref16}. The linearized equation is 


(1 + ecos(v))d@" 
3 
+ 5 (2 cos(2wv) ô 


sin(2wv) e) 4 


(23) 


w +1 


Applying the Floquet theory find the monodromy 
matrix (A) for this system, one is obtained numerically by 
integrating a periodic orbit (2m) and the variation 
equations (Eq. (23)). 


A good numerical method to refine the closure of the 
orbit is essential to obtain an accurate monodromy matrix 
(A), which is obtained to analyze to stability with respect 
to small perturbations (ôq) of the orientation angle @. It is 
possible to analyze the stability using the Floquet theory, 
since this equation is a differential equation of the second 
order. The linearized equation constricts this analysis to 
small variations of ọ. The stability monodromy matrix (A) 
has some important properties, which are ({13] and [19]): 


1. det(A)=1; 
2. {A, = à; à, = A“*}eigenvalues; 


The indicator of stability 2—|Tr(A)| is given by 
Tr(A) between 0 and 2 (0<2-|Tr(A)| S 2) 
\cite{ref13}, where Tr(A) means Trace of the matrix A. 
Positive values between the described boundaries 
correspond to stable solutions (linear approximation), 
negative values correspond to instability and zero values 
correspond to critical cases. The conditions of stability for 
œw are demonstrated in Figures 8 -- 18, where the positive 
values of w correspond to direct rotations (direction of the 
orbital motion). 


Table l - Intervals of stability. 


[0, 0.2956] U [0.9063, 0.9334] 


(0.7821, 0.8333] 


[0, 0.0966] U [0.6824, 0.7690] 


(0.6328, 0.7542] 


-3 [0, 0.7895] 
[0, 0.6370] 


275 ee 
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-2.25 [0, 0.1935] 


3 
e 
c 


[0, 0.9999] 


O E a 
0.5 No Solutions found 
A 


[0.0, 0.1338] 


[0, 0.5539] U [0.9778, 0.9789] 
0, 0.4207] 


The regions of stability with solutions (Figures 8-18 





and Table 1) have been shown for a number of w E 
[—4,4]. When w = —1 there is no solution and when w = 
1 there is a small stable range (Table 1). The Table 1 
shows the complete set of stable solutions for the 
monodromy matrix, as a function of the eccentricity and of 
the true anomaly. 


www.ijaers.com 











Fig.9: Stability curve for w = {-3; -=;-2;- = 





Fig.10: Stability curve for w = {— =| 


Fig.11: Stability curve for w = fo; =} 


4 
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Fig.12: Stability curve for w = fi; ž, n , Pi 
p 2 2 Fig.16: Stability curve for w = {—4; —3; —2}. 
2-]|Tr Al 
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Fig.13: Stability curve for w = {2; =; >; =| 
ae: 3 Hi. 
Fig.18: Stability curve for w = { 0; 1}. 
Fig.14: Stability curve for w = {3; = 57; =; 4} 


The stability analysis suggests the system's viability. 
Some parameters for the system are shown in monodromy 


matrix. 





Hl. FORCES IN TETHERS 


The force on the tethers can be calculated using the 
relation of the forces involved in the problem, applied to 


the tether, 
ws mı — 4 
Mı = —Ho FE ri +T (24) 
1 


3 The force on the tether for a system composed of: 
mı =m, = 1000 kg, lọ = 100km, p = 7000km, the 


Fig.15: Stability curve for w = { aC 
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direction is tether length (/) and sense m, to m3, and Tis T 
magnitude of the tether force, solving the Eq.(16) with 
respect to l, it is possible to show that the cable suffers a 
large variation in tension and that these values are periodic o l 
and, in some cases, increase with grows to the eccentricity = 
(Figures 19 -- 32). 


w=0 and e=0 





+— vfrad] 


wU =s 


0.30 l l 

Fig.22: Tether Force T for w = 1. 
0.25 
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i] T m IE $ pe 
i a = am 


Fig.19: Tether Force T for w = 0 ande = 0. 
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ant Fig.23: Tether Force T for w = 2. 
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Fig.20: Tether Force T for w = {= 1} and e=0. 
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w=0 Fig.24: Tether Force T for w = 3. 
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Fig.21: Tether Force T for w = Q. Fig.25: Tether Force T for w = = 
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Fig.26: Tether Force T for w = = 
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Fig.27: Tether Force T for w = -< 
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Fig.28: Tether Force T for w = — = 


w=0 


Fig.29: Tether Force (T) with w = 0 for several 


eccentricities (e). 
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Fig.30: Tether Force (T) with w = 1 for several 


eccentricities (e). 


100097 


im 


Fig.31: Tether Force (T) with w = - for several 


eccentricities (e). 





Fig.32: Tether Force (T) with w = —4 for several 


eccentricities (e). 
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IV. CONCLUSION 


The uniform rotations of a dumbbell and with several 


possibilities are considered in the present study, as well as 


the stability analysis and the viable control laws. In some 


cases are possible obtain solutions closes-form, in other 


cases only numerical solutions for the control of the tether 


systems are available. The necessary conditions of stability 


for uniform rotations were analyzed using the parameters 


œw and the eccentricity of the orbit, generating the control 


laws for the cable length. 
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